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We expand our pervious numerical study of the properties of the stellar velocity distribution within the disk of a two- 
armed spiral galaxy by considering spiral stellar density waves with different global Fourier amplitudes, C2 . We confirm 
our previous conclusion that the ratio (t\ : CJ2 of smallest versus largest principal axes of the stellar velocity ellipsoid 
becomes abnormally small near the outer edges of the stellar spiral arms. The extent to which the stellar velocity ellipsoid 
is elongated (as compared to the unperturbed value typical for the axisymmetric disk) increases with the strength of the 
spiral density wave. In particular, the C2 ~ 0.06 spiral can decrease the unperturbed value of cri ; (J2 by 20%, while the 
C2 ~ 0.13 spiral can decrease the unperturbed cri : CT2 by a factor of 3. The abnormally small values of the ai : a2 ratio 
can potentially be used to track the position of stellar spiral density waves. The a^^ : a^r ratio is characterized by a more 
complex behaviour and exhibits less definite minima near the outer edges of the spiral arms. We find that the epicycle 
approximation is violated near the spiral arms and cannot be used in spiral galaxies with C2 <; 0.05 — 0.06 or in galaxies 
with the amplitude of the spiral stellar density wave (relative to the unperturbed background) of order 0.1 or greater. 



© 0000 WILEY- VCH Vcrlag GmhH & Co. KGaA, Wcinhcim 



1 Introduction 

It is well known that the local velocities of stars in galac- 
tic disks are non-isotropic. The velocity distribution of stars 
is usually described by the Schwarzschild distribution func- 
tion, which is a usual Gaussian distribution but with differ- 
ent squared velocity dispersions cr^^, (7|^, and in the 
radial (r), azimuthal (</>), and vertical (z) coordinate direc- 
tions. In addition, the cross-correlated velocity dispersions 
(e.g. cr^^ ) may be non-zero and the principal axes of the 



velocity dispersion tensor (i, j 
aligned with the coordinate axes. 



r, 0, z) may not be 



The gravitational field of stellar spiral arms can also 

give s rise to anomalies in the stellar veloci ty distribution 

(e.g.lMavoiiri970l:lMinchev & Ouillenll2007h . Recentlv. lVorobvov & Thels^ 
( 2008h have studied numerically the dynamics of stellar disks| 
with a saturated two-armed spiral structure using the Boltz- 
mann moment equation up to second order. They find that 
the spiral gravitational field introduces large-scale non-circulaij 
motions in stellar orbits near the outer edges of the spiral 
arms. These stellar streams are responsible for the peculiar 
properties of the stellar velocity distribution such as non- 
zero vertex deviations, abnormally low values of the ratio 
a\ : (72 of smallest versus largest principal axes of the stellar 



In an axysimmeti-ic, stationary disk stellar orbits ai-e nearly|velocity dispersion tensor, and lai-ge deviations of the ratio 



circular and characterized by epicycle motions with small 
radial amplitudes with respect to the guiding center. In such 
a disk the Oort ratio = cr'^^/crf.^ can be expressed as 



= —B/ {A—B), where A and B are the usual Oort con- 
stants. For the most interesting case of a flat rotation curve, 
this gives = 0.5. However, the Solar Neighbourhood 
value for the Oort ratio has been demonstrated by almost 
all ob servers to be less than 0.5 (see e.g. Kerr & Lynden-BelH 
19861) . A possible explanation for this discrepancy between 
X^p and is the non-axysimmetric gravitational perturba- 
tion from the Galactic bar that causes significant perturba- 
tions to the velocity moments (e.g. iMiihlbauer & Dehnenl 
2003h . In addition, the Oort con stants can be affected by the 
bar's potential (e.g. lMinchev et al.. 2007). 



""00 • '^rr from those predicted by the epicycle approxima- 
tion {Xl^). 

In this paper, we expand our previous numerical study 
of the stellar velocity distribution in spiral galaxies with a 
two-armed structure ( Vorobyov & T heis. 2008|) by consid- 
ering stellar spiral density waves of different amplitude. The 
paper is organized as follows. The description of our model 
stellar disk is given in Section |2] The peculiar properties of 
the stellar velocity ellipsoid and the violation of the epicy- 
cle approximation are discussed in Sections[3]and|4] respec- 
tively. The main conclusions are summarized in Section|5] 
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2 Model description 

We use the BEADS-2D co de that is designed to stud y the 
dynamics of stellar disks ( Vorobyov & Theisi 20061) . The 
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BEADS-2D code is a finite-difference numerical code that 
solves the Boltzmann moment equations up to second order 
in the thin-disc approximation on the polar grid (r, </)). More ^ 
specifically, the BEADS-2D code solves for the stellar sur- § 
face density S, mean radial and azimuthal stellar velocities g 
Ur and u^, and stellar velocity dispersion tensor Oij. The ro 
latter includes the squared radial and azimuthal stellar ve- "5 
locity dispersions a^^ and cr"^^, respectively, and the mixed 
velocity dispersion a^^ . We close the system of Boltzmann 
moment equations by adopting the zero-heat-flux approxi- 
mation. For the main equations and necessary tests we refer 
the reader to lVorobyov & TheisI (I2006ll2008h . 

Our model galaxy consists of a thin, self-gravitating stel-| 
lar disc embedded in a static dark matter halo. The initial 
surface density of stars is axisymmetric and is distributed 
exponentially according to 

E(r) = Eoexp(-r/rd), (1) 

with a radial scale length of 4 kpc and central surface 
density Sq = 10"^ Mq pc^^. The gravitational potential of 
the stellar disk is calculated by solvin g for the Poisson in- 
tegral using the convolution theorem dBinnev & Tremaine , 
19871 section 2.8). 

The initial mean rotational (azimuthal) velocity of stars 
in the disc is chosen according to 









c^^o 13 '^^^^^^^^ 



- 

V'"flat 

The transition radius between an inner region of rigid ro- 
tation and a flat rotation in the outer part is given by rfiat, 
which we set to 3 kpc. The smoothness of the transition is 
controlled by the parameter riy, set to 3. The velocity at in- 
finity (Uoo) is set to 208 km s^^. 

The radial component of the stellar velocity dispersion 
is obtained from the relation (!„ = 3.36(5sGS/k for a 
given value of the Toomre parameter Qs- Here, k is the 
epicycle frequency. We assume that throughout most of the 
disc Qs is constant and equal 1 .3 but is steeply increasing 
with radius at r > 25 kpc. The azimuthal component of 
the velocity dispersion ct^^ is determined from X'^^ using 
the initial rotation velocity of stars in the Oort constants. 
Once the rotation curve and the radial profiles of the stellar 
surface density and velocity dispersions are fixed, the dark 
matter halo potential can be derived from the steady-state 
momentum equation for t he azimuthal velocity of stars (see 
Vorobvov & TheisLl2008l for details). 

In the current simulations, the numerical resolution has 
512x512 grid zones that are equally spaced in the (/)-direction| 
and logarithmically spaced in the r-direction. The inner and 
outer reflecting boundaries are at rin = 0.2 kpc and rout = 
35 kpc, respectively. The total disk mass Mst = 10^^ Mq. 
We extend the outer computational boundary far enough 
to ensure a good radial resolution in the region of interest 
within the inner 20 kpc. For instance, the radial resolution 
at 1 kpc and 20 kpc is approximately 10 pc and 150 pc, re- 
spectively. 
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Fig. 1 Relative stellar density perturbation (black contour 
lines) superimposed on the grey-scale image of the cti : a2 
distribution at t = 1.5 Gyr (left) and t = 1.6 Gyr (right). 
The corresponding global Fourier amplitudes of the m=2 
mode (C2) are indicated in each panel. 

3 Peculiar shape of the stellar velocity 
ellipsoid 

In the beginning of our simulations we add a small (< 10^^) 
random perturbation to the initially axisymmetric surface 
density distribution of stars. Since the model stellar disc is 
characterized by Qs = 1-3, it becomes vigorously unsta- 
ble to the growth of non-axisymmetric gravitational insta- 
bilities. The initial parameters of our model are chosen to 
favour the growth of a two-armed spiral pattern. The most 
likely physical interpretation for the growth o f a spiral struc- 
ture in our model disk is swing amplification (Voro byov & Theisl j 
I2OO6 , 200 8). The black contour lines in Fig.[T]show the rel- 
ative stellar density perturbation defined as 



C(r-,0) 



Si„(r) 



(3) 



where Sin ('') is the initial axisymmetric stellar density dis- 
tribution in the disk and E(r, (f) is the current stellar density 
distribution. Only positive density perturbations are plotted. 
The left panel in Fig.[T|corresponds tot = 1.5 Gyr since the 
beginning of our simulations and the minimum and maxi- 
mum contour levels correspond to the relative perturbations 
of 0.03 and 0.23, respectively. The right panel shows the 
relative stellar density perturbations at t — 1.6 Gyr The 
minimum and maximum contour levels are C = 0.05 and 
C = 0.75, respectively. To quantify the strength of the spiral 
arms, we use the global Fourier amplitudes defined as 

2tt 



Cn.it) = 
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The dominant m — 2 mode in the left/right panels of Fig. [T] 
is characterized by the C2 = 0.06 and C2 = 0.13, respec- 
tively. 

The velocity dispersion tensor af^ is often non-diagonal 
in the local coordinate system (r, 0, z). The principal axes 
of a diagonalized velocity dispersion tensor form an imagi- 
nary ellipsoidal surface that is called the velocity ellipsoid. 
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Fig. 2 Radial profiles of the cti : (72 ratio (solid lines), 
o'<p<p ■ cTrr ratio (dashed lines), and relative stellar density 
perturbation Q (shaded area) obtained by taking a narrow 
radial cut at = 145° in the top panel of Fig. [T] (top) and 
(j) = 110° in the bottom panel of Fig.[T|(bottom). The dotted 
lines present the unperturbed initial values of ai : a2- The 
corresponding strength of the m = 2 mode is indicated in 
each panel. 



To better illustrate this phenomenon, we take a narrow 
radial cut in Fig.[I]at = 145° (top) and (f) = 110° (bot- 
tom), and plot the resulting radial distributions of ( (shaded 
area) and ai : (72 (solid lines) in Fig. |2l The angle (j> is 
counted counterclockwise from the positive horizontal axis 
of Fig. [U The dotted line shows the initial unperturbed val- 
ues of (7i : (72. It is evident that the weaker C2 = 0.06 spiral 
(top) causes only a mild perturbation to the initial shape of 
the stellar velocity ellipsoids. However, there is a notice- 
able decrease in the values of cti : (72 near the outer edge 
of the innermost (and strongest) spiral arm, indicating that 
the shape of the stellar velocity ellipsoid becomes abnor- 
mally elongated there. As the global Fourier amplitude of 
the m — 2 mode increases to C2 = 0.13 (bottom), a con- 
siderable change in the initial shape of the stellar velocity 
ellipsoid becomes evident. Deep minima in the radial distri- 
bution of (7i : (72 have occurred near the outer edges of the 
spiral arms. We conclude that the peculiar elongated shape 
of the stellar velocity ellipsoids can be used as a tracer of the 
stellar spiral density wave if the latter is sufficiently strong, 
C2 ^ 0.1. 

It may be tempting to use the ratio a^^ : arr instead 
of (7i : (72 in order to track the position of stellar spiral 
density waves. However, we find that the latter ratio may 
not be as good a diagnostic tool as the former ratio. Indeed, 
the dashed lines in Fig. |2] show the radial distribution of 
the CT0^ : CTrr ratio taken along the same azimuthal angles 
as in the case of the ai : 02 ratio. It is evident that the 
1700 : <7rr ratio also exhibits local minima near the outer 
edges of the spiral arms but these depressions may not be as 
deep as those of the a\ : (J2 ratio. 



The available measurements in the solar vicinity indicate 
that the ratio a\ : G2 of smallest versus largest principal 
axes of stellar velocity ellipsoids in the disk plane does not 
vary significantly with the B — V co lour or with the age 
stellar populations ( Dehnen & Binneviil998h . On the other 
hand, the : (J2 ratio is expected to have deep minima near 
the outer edges of spiral arms and this peculiar, elongated 
shape of the stellar velocity ellipsoid ca n be used to track the 
positi on of spiral stellar density waves d Vorobvov & Theis , 
20081). In the present paper we corroborate our previous re- 
sults by considering spiral stellar density waves with dif- 
ferent global Fourier amplitudes. The grey-scale images in 
Fig. [T] show the spatial distribution of a\ : 02 obtained at 
t Gyr (left), when the global Fourier amplitude of the 

m — 2 mode is C2 = 0.06, and at t = 1.6 Gyr (right), when 
C2 has reached a maximum value of 0.13. It is seen that the 
stronger spiral produces a larger response in the shape of the 
stellar velocity ellipsoids. For instance, the minimum values 
of(7i : (72 in the C2 = 0.06 disk are about 0.5 — 0.6, but they 
decrease to 0.2 — 0.3 in the C2 = 0.13 disk. In both cases 
the shape of the stellar velocity ellipsoids is most elongated 
near the outer edges of the spiral arms, with the exception 
of a small region near the central bar. 



4 Violation of the epicycle approximation 

It is often observationally difficult to measure both radial 
and azimuthal velocity dispersions of stars. Hence, it is tempt-| 
ing to use the Oort ratio to derive either of the two disper- 
sions from the known one and the circular speed of stars 

-B \ f r du^ 
Un dr 



1 



(5) 



A-B 2 

One has to keep in mind that equation (|5]l may be valid only 
in the epicycle approximation (e.g Kuijken & Tremainei , ll994l) [ 

To test the validity of equation Q, we calculate the fol- 
lowing quantity that measured the relative deviation from 
the epicycle approximation 

i{r,<f>) = — , (6) 

^(p(f) • (7 rr 

where the ratio (7^0 : cth is calculated from the model's 
known velocity dispersions and is calculated using equa-| 
tion ^ and u^, as a proxy for Uc- The resulting radial pro- 
files of ^ taken along the azimuthal angles (j) ~ 145° (top) 
and (j) — 110° (bottom) are plotted in Fig. [3] by the solid 
lines. The shaded area shows the positive stellar density per- 
turbation C for the C2 = 0.06 spiral (top) and C2 = 0.13 
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Fig. 3 Radial profiles of ^ (solid lines) and relative stellar 
density perturbation C, (shaded area) obtained by taking a 
narrow radial cut at (/> = 145° in the top panel of Fig. [T] 
(top) and = 1 10° in the bottom panel of Fig. [T](bottom). 
The corresponding strength of the m = 2 mode is indicated 
in each panel. 

spiral (bottom), respectively. It is seen that the epicycle ap- 
proximation is violated near the outer edges of the spiral 
arms. The extent to which the epicycle approximation is vi- 
olated depends on the strength of the spiral structure. For in- 
stance, the C2 = 0.06 spiral shows a moderate deviation of 
order 50% near the position of the innermost spiral arm and 
a mild deviation of order 20% near the second spiral arm. 
As the strength of the dominant m = 2 mode has increased 
to C2 = 0.13, the deviation from the epicycle approxima- 
tion begins to exceed 100%. It is therefore very risky to use 
the epicycle approximation to derive velocity dispersions, 
when the global Fourier amplitude of the dominant spiral 
mode exceeds 0.05 — 0.06. Figure[3]indicates that the same 
criterion can be expressed in terms of ^ - the relative am- 
plitude of the spiral stellar density wave should not exceed 
0.1. 



- The stellar spiral density wave produces peculiar signa- 
tures in the velocity distribution of the underlying stellar 
disk. For instance, the ratio ai : (T2 of the smallest ver- 
sus largest principal axes of the stellar velocity ellipsoid 
becomes abnormally small (as compared to the unper- 
turbed value typical for an axysimmetric stellar disk) 
near the outer edges of the spiral arms. The degree to 
which the stellar velocity ellipsoid becomes elongated 
depends on the strength of the spiral arms (as defined 
by the global Fourier amplitudes. Cm). For instance, the 
C2 = 0.06 spiral can reduce the initial (unperturbed) 
value of CTi : 0-2 from 0.75 to 0.6, while the C2 — 0.13 
spiral can reduce the initial cti : (T2 by a factor of three. 
The abnormally small values of the di : 0-2 ratio can 
potentially be used to track the position of stellar spiral 
density waves. 

- The epicycle approximation is violated near the spiral 
arms of galaxies with the global Fourier amplitude of 
the dominant m — 2 mode of order C2 — 0.05 — 0.06 
or greater. This corresponds to the stellar density pertur- 
bation (relative to the unperturbed background disk) of 
order 0.1 or greater 

- The (700 : (Trr ratio of the azimuthal to radial velocity 
dispersions shows less definite minima near the outer 
edges of the spiral arms. 
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5 Conclusions 

We have studied numerically the stellar velocity distribution 
in a two-armed spiral galaxy with a varying strength of the 
dominant m = 2 mode. The development of the spiral struc- 
ture and its time evolution are followed self-consistently by 
solving the Boltzmann moment equations up to second or- 
der in the thin-disk approximation. We find the following: 
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